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Abstract 

We calculate the chiral and thermal susceptibilities for two confining Dyson- 
Schwinger equation models of QCD with two light flavours, a quantitative 
analysis of which yields the critical exponents, (3 and 5, that characterise the 
second-order chiral symmetry restoration transition. The method itself is of 
interest, minimising the influence of numerical noise in the calculation of the 
order parameter for chiral symmetry breaking near the critical temperature. 
For the more realistic of the two models we find: T c ~ 153 MeV, and the non- 
mean-field values: (3 = 0.46 ± 0.04, 5 = 4.3 ± 0.3 and 1/(J38) = 0.54 ± 0.05, 
which we discuss in comparison with the results of other models. 
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I. INTRODUCTION 



Phase transitions are characterised by the behaviour of an order parameter, (X), the 
expectation value of an operator. In the ordered phase of a system: (X) ^ 0, whereas in 
the disordered phase (X) = 0. A phase transition is first-order if (X) — > discontinuously, 
whereas it is second-order if (X) — > continuously. For a second-order transition, the length- 
scale associated with correlations in the system diverges as (X) — > and one can define a 
range of critical exponents that characterise the behaviour of certain macroscopic properties 
at the transition point. For example, in a system that is ferromagnetic for temperatures less 
than some critical value, T c , the magnetisation, M, in the absence of an external magnetic 
field, behaves as M oc {T c — Tf for T ~ T~, where (5 is the critical exponent. At the critical 
temperature the behaviour of the magnetisation in the presence of an external field, h — > + , 
defines another critical exponent, 6: M oc h^' s \ In a system that can be described by mean 
field theory these critical exponents are 

/5 MF = 0.5, 5 MF = 3.0. (1) 

Equilibrium, second-order phase transitions can be analysed using the renormalisation 
group, which leads to scaling laws that reduce the number of independent critical exponents 
to just two: j3 and 5 [|TJ. It is widely conjectured that the values of these exponents are 
fully determined by the dimension of space and the nature of the order parameter. This 
is the notion of universality; i.e., that the critical exponents are independent of a theory's 
microscopic details and hence all theories can be grouped into a much smaller number of 
universality classes according to the values of their critical exponents. If this is the case, the 
behaviour of a complicated theory near criticality is completely determined by the behaviour 
of a simpler theory in the same universality class. So, when presented with an apparently 
complicated theory, the problem is reduced to only that of establishing its universality class. 

Quantum chromodynamics is an asymptotically free theory; i.e., there is an intrinsic, 
renormalisation-induced mass-scale, A QCD , and for squared momentum transfer Q 2 ^> Aqcd, 
the interactions between quarks and gluons are weaker than Coulombic: as(Q 2 ) — ► as 
Q 2 — > oo. The study of QCD at finite temperature and baryon number density proceeds via 
the introduction of the intensive variables: temperature, T; and quark chemical potential, \i. 
These are additional mass-scales, with which the coupling can run and hence, for T 3> Aqcd 
and/or \i ^> A QCD , as(Q 2 = 0, T, //) ~ 0. It follows that, at finite temperature and/or baryon 
number density, there is a phase of QCD in which quarks and gluons are weakly interacting, 
irrespective of the momentum transfer ||; i.e., a quark-gluon plasma phase. Such a phase 
of matter existed approximately one microsecond after the big-bang. 

At T, jj, — the strong interaction is characterised by confinement and dynamical chiral 
symmetry breaking (DCSB), effects which are tied to the behaviour of as(Q 2 ) at small-Q 2 ; 
i.e., its long-range behaviour. In a phase of QCD in which the coupling is uniformly small 
for all Q 2 , these effects are absent and the nature of the strong interaction spectrum is 
qualitatively different. 

The path followed in the transition to the plasma is also important because it deter- 
mines some observational consequences of the plasma's existence. For example [[|, the 
time-scale for the expansion of the early universe: ~ 10~ 5 s, is large compared with the 
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natural time-scale in QCD: 1/Aqcd ~ lfm/c ~ 10 _23 s, hence thermal equilibrium is main- 
tained throughout the QCD transition. Therefore if the transition is second-order the ratio 
B := baryon- number/entropy, remains unchanged from that value attained at an earlier 
stage in the universe's evolution. However, a first-order transition would be accompanied by 
a large increase in entropy density and therefore a reduction in B after the transition. Hence 
the order of the QCD transition constrains the mechanism for baryon number generation in 
models describing the formation of the universe, since with a second-order transition this 
mechanism is only required to produce the presently observed value of B and need not al- 
low for dilution. In the absence of quarks, QCD has a first-order deconfinement transition, 
while with three or four massless quarks a first-order chiral symmetry restoration transition 
is expected |§. What of the realistic case with two light quark flavours? 

Based on the global chiral symmetry of QCD with two light quark flavours, it has been 
argued that this theory and the N = 4 Heisenberg magnet are in the same universality 
class. As a field theory, the N = 4 Heisenberg magnet is characterised by an interaction of 
the form 

E{|A 2 (x) + iA 4 0^)}, (2) 
i=i 

where fi 2 is a function of temperature: /i 2 > at or above the critical temperature, Tf , 
but [j 2 < for T < T^ 1 . If the interaction strength, A, depends smoothly on T and remains 
positive then, for T < , the classical minimum of this potential is at 

^ = ^!> . (3) 

This model is familiar as the nonlinear cr-model, often used to describe low-energy phenom- 
ena in QCD. It has been explored thoroughly and has a second order phase transition with 
critical exponents ||J] 

(3 H = 0.38 ±0.01, 5 H = 4.82 ±0.05. (4) 

One can examine the hypothesis that this model and QCD with two light quark flavours 
are in the same universality class via numerical simulations. Such studies on an 8 3 x 4 lattice 
suggest a second-order chiral phase transition with critical exponents |5[] 

/3 lat = 0.30 ± 0.08 , 5 lat = 4.3 ± 0.5 (5) 

but do not decide the question.^] These results were obtained through an analysis of the 
chiral and thermal susceptibilities; a technique that can be applied in the study of any 



1 A review || of results from more recent simulations on larger lattices with lighter quarks reports 
a significant dependence of these critical exponents on the lattice volume but with their product 
approximately constant. A value of 5 ~ 1 is obtained, which is characteristic of a first-order 
transition. These unexpected results might be artefacts of finite lattice spacing because introducing 
light dynamical quarks drives the simulations to stronger coupling and hence coarser lattices. 
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theory. Herein we illustrate the method via an analysis of two Dyson- Schwinger equation 
(DSE) models of QCD, which also allows us to explore the hypothesis further. 

Dyson-Schwinger equations provide a renormalisable, nonperturbative, continuum frame- 
work for the exploration of strong interaction effects and have been used extensively at 
T = 0] in the study of confinement and DCSB, and in the calculation of hadron observ- 
ables pHlOf . They have recently fTl"l , |T2l found successful application at T 7^ and it is these 
two models that we employ as exemplars herein. In Sec. II we describe the models and in 
Sec. Ill the analysis of their chiral and thermal susceptibilities, and the evaluation of the 
associated critical exponents. We summarise and conclude in Sec. IV. 



II. TWO MODELS 

Using a Euclidean metric, with {7^,7^} = 28^ v and 7^ = 7^, the renormalised dressed- 
quark propagator at T 7^ takes the form 

S(p Uk ) ■= -%1 ■ pvA(Puj k ) - il4U k <j c (Pu k ) + cr B (Pw k ) , (6) 

where (p Uk ) '■= (p, oj k ) with uj k = (2k + 1) nT the fermion Matsubara frequency, and a^(p iUk ), 
T = A,B, C, are functions only of \p\ 2 and The propagator is obtained as a solution of 
the quark DSE 

S' 1 (Pu lk ) := ij ■ p A(p UJk ) +ry4^fcCG%J + B(p L , k ) (7) 
= Zztf ■ p + Z 2 (i'jiUJk + m bm ) + £'(p„J, (8) 

mb m is the Lagrangian current-quark bare mass and the regularised self energy is 

wit) — ^7 ' P^A(P^k) + ^74 ^c(P^k) + , (9) 

with 

/•A 

^(P«J = / i9 2 D^(p- q,uj k -uJi)\ti[P^ ll S(q UJt )T u (q UJl ;p L , k )} , (10) 

J l,q 

where: T = A,B,C\ V A := -(Zf_/\p\ 2 )ii-p, V B := Z u V c := -(Zjoo^; and 
J [q := T Y^"iZ-oo f A d 3 q/(2ir) 3 , with J A a mnemonic to represent a translationally invariant 
regular isat ion of the integral and A the regularisation mass-scale. In Eq. flip]), T u (q UJl ;p UJk ) 
is the renormalised dressed-quark-gluon vertex and D^(p,Qk) is the renormalised dressed- 
gluon propagator. (Q k = 2k nT is the boson Matsubara frequency.) 
In renormalising the quark DSE we require that 

S ~HPuo)\ lpl2+UJ 2 =e =il-p + ilAU Q + m R , (11) 

which entails that the renormalisation constants are 

z 2 A (c,A) = i-£' A (c ;A), (12) 

Z 2 (C,A) = 1-£' C (C ;A), (13) 
m R (C) = Z 2 m hra (K 2 ) + E' fl (C ; A), (14) 
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where (C, ) 2 := C 2 — an d the renormalised self energies are 

Hp^O = ^ + S^ fc ;A)-E^(C ;A), (15) 

T = A, B, C, U = 1 = £c and £ B = m R {(). 

So far no approximations or truncations have been made but to continue we must know 
the form of T„(q Ul ;p Uk ) and D^ u (p, Q k ) in Eq. fllpp. These Schwinger functions satisfy DSEs. 
However, the study of those equations is rudimentary even at T = and there are no studies 
for T ^ 0. To proceed we use the T = results as a qualitative guide and employ exploratory 
Ansdtze for r t/ (g aJi ;p a;fc ) and D^^flk). This is where model parameters enter. 



The structure of the dressed fermion-gauge-boson vertex has been much considered |L3 
As a connected, irreducible three-point function it should be free of light-cone singularities 
in covariant gauges; i.e., it should be regular at (p — q) 2 + {oJk — ^i) 2 = 0. A range of Ansdtze 
with this property have been proposed and employed (TTJJ and it has become clear that the 
judicious use of the rainbow truncation 

r u (q Ul ;p Wh ) = 7„ (16) 



in Landau gauge provides reliable results [15]. This is the Ansatz employed in Refs. [11,12 



and we use it herein. With this truncation a mutually consistent constraint is Z\ = Z2 and 

With r j/ (g aJ; ; p Wfc ) regular, the analytic properties of the kernel in the quark DSE are 
determined by those of D^ v (pci k ), which in Landau gauge has the general form 

g 2 D^(pn k ) = P^(pn k )A F (p nk ) + Pj> n JA G (p n J , (17) 



0; fi and/or v 



CteiJ = U,-?pw = ij = 1,2,3, < 18 > 

with P^ u (pn k ) + P^uiPnJ = <W ~ mVELi PaPa] V, v = 1, . . . , 4. A "Debye-mass" for 
the gluon appears as a T-dependent contribution to A F . Considering D^ u (k) at T = 0, 
a perturbative analysis at two-loop order provides a quantitatively reliable estimate for 
k 2 > 1-2 GeV 2 , with higher order terms providing corrections of only ~ 10%. However, 
for k 2 < 1 GeV 2 nonperturbative methods are necessary. Studies of the gluon DSE in 



axial gauge [pl|, where ghost contributions are absent, or in Landau gauge [ITU, when their 
contributions are small, indicate that D liv {k) is significantly enhanced in the vicinity of 
k 2 = relative to a free gauge-boson propagator, and that the enhancement persists to 
k 2 ~ lGeV 2 . Due to the truncations involved these studies are not quantitatively reliable 
but this behaviour has been modelled successfully as a distribution located in the vicinity 
of k 2 = il5|JH. 



A. Infrared- dominant Model 

A particularly simple and illustratively useful model is obtained with 
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A F (p nk ) = A G ( Pnk ) = 2tt 3 ^S k0 6 3 (p) , (19) 



which is a generalisation to T 7^ of the model introduced in Ref. [[HJ], where rj m 1.06 GeV 
is a mass-scale parameter fixed by fitting 7T- and p-meson masses. As an infrared-dominant 
model Eq. (|T^) does not represent well the behaviour D^ u (pa k ) away fromp 2 ^ ~ 0, and hence 
there are some model-dependent artefacts. However, these artefacts are easily identified and, 
because of its simplicity, the model has provided a useful means of elucidating many of the 
qualitative features of more sophisticated Ansdtze. 

Using Eqs. (|l~6l) and (|i~9"D the quark DSE is ultraviolet-finite, the cutoff can be removed 
and the renormalisation point taken to infinity, so that Eq. (||) becomes the algebraic equa- 
tions 

v 2 m 2 = B 4 + mB 3 + Upl k - - m 2 ) B 2 - m (2 rf + m 2 + 4p 2 J B , (20) 



with Z2 = 1 = Z2 and m = = mb m : m = defines the chiral limit. This DSE- model 
of QCD has coincident, second-order deconfinement and chiral symmetry restoration phase 
transitions at a critical temperature T C IR 0.16 1] flT2 . 



B. Ultraviolet-improved Model 

An improvement over Eq. (|l9|) is obtained by correcting the \aige-pf lk behaviour so as to 
better represent the interaction at short-distances. The one-parameter model 



A F (pn k ) = V(p nk ;m D ) 
A G (pn h ) = V(p nk ;0), 



V(p nk ;m) :=f vr 2 



2vr 1 _ e h(|pl 2 +^+m 2 )/(4m2 

— m 2 t 5 0k 5 3 (p) + — 



(22) 
(23) 

(24) 



where m 2 ^ 



'8/3)7r 2 T 2 is the perturbatively evaluated "Debye-mass"0, achieves this. This 
gluon propagator provides a generalisation to T 7^ of the model explored in Ref. |18j where 
the parameter m t is a mass-scale that marks the boundary between the perturbative and 
nonperturbative domains. The value m t = 0.69 GeV = 1/0.29 fm is fixed by requiring a 
good description of a range of tt- and p-meson properties. In this case the DSE yields a pair 
of coupled, nonlinear integral equations that must be solved subject to the renormalisation 
boundary conditions, and m# = defines the chiral limit. This model also has coincident, 
second-order deconfinement and chiral symmetry restoration transitions, with the critical 
temperature T c uv 



0.15GeV |11]. 



2 The influence of the Debye-mass on finite-T observables is qualitatively unimportant, even in the 
vicinity of the chiral symmetry restoration transition. The ratio of the coefficients in the two terms 
in Eq. ([24]) is such that the long-range effects associated with §ok^ 3 (p) are completely cancelled at 
short-distances; i.e., for |x| 2 m 2 <C 1. 
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III. CHIRAL AND THERMAL SUSCEPTIBILITIES 



In the study of dynamical chiral symmetry breaking an order parameter often used is 
the quark condensate, (qq)^- In QCD in the chiral limit this order parameter is defined via 
the quark propagator fl~5H : 

- (qq) c := N c Jim Z 4 ((, A) WT^T W^Wl v > ( 25 ) 

for each massless quark flavour, where the subscript "0" denotes that the scalar functions: 
A , Bq, C , are obtained as solutions of Eq. (||) in the chiral limit, and Z 4 (C,A) is the 
mass renormalisation constant: Z 4 (£,A)?71r(C) = Z 2 ((, A) mb m (A). The functions have an 
implicit ^-dependence. From Eq. ( p5| ) it is clear that an equivalent order parameter for the 
chiral transition is 



X := B (p = 0,u ) 



(26) 



which was used in Refs. 11,12 



Thus the zeroth Matsubara mode determines the character 
of the chiral phase transition, a conjecture explored in Ref. 20 . 

To accurately characterise the chiral symmetry restoration transitions in the two models 
introduced above, we examine closely the chiral and thermal susceptibilities and their scaling 
behaviour near the critical point. This allows a determination of the critical temperature, 
T c , and exponents (3 and S, as we explain in the appendix. In the notation of the appendix, 
the "magnetisation" is 



M(t,h) := B(p = 0,u;o) 



(27) 



i.e., the value in the infrared of the scalar piece of the quark self energy obtained as the rriR- 
and T-dependent solution of Eq. (|8|). 



A. Critical Exponents of the Infrared-dominant Model 

In the chiral limit, Eq. (^) has the Nambu-Goldstone mode solution 



and hence 



B(pu, k ) = { ^ 4 ^ fc ' p2j 1 fc< " 4 (28) 



, otherwise 



2 



2 v 2 < r -£ 



C(P»k) = { if,, /l ; 2r,A (29) 



i(l + Jl + ^h) , otherwise, 



i i 



M(t,0) = 2vr(|- + T) 2 (|--T) 2 . (30) 



From Eq. ( |3~0| ) we read that 
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T c IR = -5- « 0.159155 r? , /? IR = - . (31) 
2,71 2 

To determine 5 we use Eq. (p0|) at T = T c to obtain 

r? 2 m 2 = M(0, /i) 4 + m M(0, /i) 3 + m 2 M(0, /i) 2 - m(3r/ 2 + m 2 ) M(0, /i) (32) 

and suppose that, for m ~ 0, M(0, /i) = am 1 ^ 5 . Consistency requires 

6 m = 3 . (33) 

That the chiral symmetry restoration transition in this model is characterised by mean 
field critical exponents is not surprising because the interaction described by Eq. ( fL9|) is a 
constant in configuration space. Mean field critical exponents are also obtained in chiral 
random matrix models of QCD p0|j21[ . 



To illustrate the evaluation of the critical temperature and exponents using the chiral 
and thermal susceptibilities we use Eqs. (0), ( |A9| ) and ( |A10| ) to obtain 

2 M(T, h) T fe_ [1 — M(T, h) bj - M(T, ft) T b + 
Xh{ ,n) 2M{T,h)b_b + [l-M{T,h)b_]-b_b + -M{T,h)b_ , { ' 

87T 2 TM{T,h) 2 b 2 _-2M{T,h)b_b + h[l- M{T,h)b-]- M{T,h)hb 2 + 
XT ^ ' ' ~ 2 M(T, fc) 6_ b + [1 - M(T, fc) 6_] - 6_ 6+ - M(T, h) 6_ ' ^ ^ 

where 6± := M(T,h) ± /iT. In Fig. |I| we plot the chiral susceptibility. The temperature 
dependence is typical of this quantity, with the peak increasing in height and becoming 
narrower as h — > + ; i.e., as the external source for chiral symmetry breaking is removed. To 
understand this behaviour, recall that the chiral susceptibility is the derivative of the order 
parameter with-respect-to the explicit chiral symmetry breaking mass. Denote the typical 
mass-scale associated with DCSB by M x . For h ^> M x , explicit chiral symmetry breaking 
dominates, with the order parameter X ~ h and insensitive to T, and hence Xh ~ const. 
For h ~ M x , X begins to vary with T because the origin of its magnitude changes from the 
explicit mass to the DCSB mechanism as T passes through the pseudocritical temperature, 
Tp C . This is reflected in Xh as the appearance of a peak at T£ c . For h M x , X ~ h until 
very near T£ c when the scale of DCSB overwhelms h and X ~ M x . The change in X is rapid 
leading to the behaviour observed in Xh- The thermal susceptibility is plotted in Fig. ^| and 
has qualitatively similar features. 

In Table | we present the pseudocritical points and peak heights obtained for h in the 
scaling window, defined as the domain of h for which 

t h c 

-j— = const. ; (36) 



i.e., the values of h for which Eqs. ( A15 ) and ( A16 ) are valid. Based on Eqs. ( |A17| ) and 



( |A18|) , using the tabulated values, one obtains z 1 ^ and zj R from linear fits to the curves: 
log Xh C ~ versus- log h and log XT- versus "l°g h, respectively. This yields 



JR _ n rrr ~IR 



0.666, z t 1K = 0.335 , (37) 



S 



and hence 0^ = 0.499 and 5* = 2.99, as listed in Table 0. These values are in excellent 
agreement with the exact (mean field) results, Eqs. ([H]) and (|33"D . With the value of 



1 - 4 R + z] R = 0.670 , (38) 



((38) 



IK 



T C IR can be obtained in a variational procedure based on Eq. (|A15|) : it is that value 
which minimises the standard deviation between log(Tp C — T C IR ) — l/(/3<5) IR log/i and a con- 
stant. This yields T C IR = 0.159155 77 again in excellent agreement with Eq. ([H]). The 
value in Table || is obtained with rj = 1.06 GeV |TJ|. Applying the same procedure to 
log(T p T c - T C IR ) - l/(pS) m hgh, yields T C IR = 0.159151 77. 



B. Critical Exponents of the Ultraviolet-improved Model 



In this case the quark DSE must be solved numerically, as in Refs. fll|jla|| . In these cal- 
culations we used a 3-momentum grid with 96 points and we renormalised at C = 9.47 GeV, 
the value at which the parameter m t (= 0.69 GeV) was fixed JTH] . The chiral and thermal 



susceptibilities for a range of values of h are plotted in Figs. | and [|, and the pseudocritical 
points and peak heights obtained for values of h in the scaling window are presented in 
Table |TTJ. 

As observed in Sec. |III A , one obtains z^ v and from linear fits to the curves logx^ - 



versus-log h and log XT C " versus "l°g K respectively. The data and fits are presented in Fig. || 
and yield 



~ uv - 0.77 ± 0.02 , 2 t uv = 0.28 ± 0.04 , (39) 



with the corresponding results for (3 and 5 listed in the first column of Table [TJ.[] For this 
model only, as a check and demonstration of consistency, the values of T c uv and l/(/35) VY 
were calculated using a variational procedure based on Eqs. ( |A15| ) and ( |A16| ): the values of 



T c uv and l/(/35) uv were varied in order to minimise the standard deviation in a linear fit to 
log(T pc — T c uv ) — 1/ (/?<5) uv log/i. The difference between using T^ c and T^ c is less than the 
error quoted in the table. 

In Ref. [JTTJ] the values of (3 and T c in the ultraviolet-improved model were calculated 
directly from the magnetisation order parameter; i.e., using Eq. (|A"6|), with the results f3 = 
0.33 ± 0.3 and T c w 152 MeV. There is a discrepancy in the value of f3. We expect that 
the result obtained herein is more accurate because our method avoids the numerical noise 
associated with establishing the precise behaviour of the order parameter in the vicinity of 
the critical temperature. 



3 Our quoted error bounds the slope of the linear fit. It is calculated from the slope of linear fits to 
the two endpoint values when they are displaced vertically, in opposite directions, by the standard 
deviation of the fit to all the tabulated results. 
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IV. SUMMARY AND CONCLUSIONS 



A primary purpose of this study was an illustration of the method by which one can calcu- 
late the critical exponents that characterise a chiral symmetry restoration transition, j3 and 8, 
using the chiral and thermal susceptibilities. For this purpose we chose two Dyson-Schwinger 

in 



equation models of two-light-flavour QCD that have been applied successfully [pl^ -pl 
phenomenological studies of QCD at finite temperature and density. The method is reliable 
and should have a wide range of application because it is more accurate in the presence 
of numerical noise than a straightforward analysis of the chiral symmetry (magnetisation) 
order parameter. 

We established that our finite temperature extension of the infrared-dominant model of 
Ref. [B5| is characterised by mean field critical exponents, listed in Table |TJ. It is therefore 
not in the universality class expected |3|,§] of two-light-flavour lattice-QCD. However, the 
critical temperature is consistent with that estimated in lattice simulations. This fits an 
emerging pattern that DSE models whose mass-scale parameters are fixed by requiring a 
good description of hadron observables at T = 0, yield a reliable estimate of the critical 
temperature for chiral symmetry restoration. It is a quantity that is not too sensitive to 
details of the model. 

Consistent with this observation, the critical temperature in the ultraviolet improved 



model of Ref. [18] also agrees with that estimated in lattice simulations. The critical ex- 



ponent S agrees with the value obtained for two-light-flavour lattice-QCD and is consistent 
with that of the N = 4 Heisenberg magnet. However, the difference between the value of 
/3 obtained in the model and that in lattice simulations is significant. It is unlikely that 
numerical errors in our study are the cause of this discrepancy. The values of the critical 
exponents, and their product, establish that this model is not mean field in character. They 
also establish that the ultraviolet-improved model, which provides a good description of low- 
energy 7r- and p-observables, is not in the same universality class as the N = 4 Heisenberg 
magnet. 

The difference between the infrared-dominant model and the ultraviolet-improved one is 
the value of m t ; i.e., the mass scale that marks the boundary between strong and weak cou- 
pling. In the limit m t — > oo, the infrared-dominant model is recovered from the ultraviolet- 
improved one: in this limit the interaction is always strong. Our results therefore demon- 
strate that the critical exponents are sensitive to the particular manner in which the theory 
makes the transition from strong to weak coupling. This should be expected since that 
evolution is a determining characteristic of the /5-function of a renormalisable theory, one 
which a chiral symmetry restoration transition must be sensitive to. 

The large-p 2 behaviour of the gluon propagator in the ultraviolet-improved model, al- 
though better than that in the infrared-dominant model, is still inadequate. Its renormal- 
isation group properties are more like those of quenched-QED than QCD because of the 
absence of the logarithmic suppression of the running coupling characteristic of asymptoti- 
cally free theories. This is corrected in the model of Ref. ]T5|] , which has more in common 
with QCD at T = and whose finite temperature properties can therefore assist in better 
understanding the details of the chiral symmetry restoration transition in two-light-flavour 
QCD. 
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APPENDIX: CRITICAL EXPONENTS FROM SUSCEPTIBILITIES 

Consider the free energy of a theory, represented by 

f = f(t,h), (Al) 

where t := T/T c — 1 is the reduced temperature and h := m/T is the explicit source of chiral 
symmetry breaking measured in units of the temperature; it is analogous to an external 
magnetic field. Since correlation lengths diverge in a second-order transition it follows that 
for t, h — > the free energy is a generalised homogeneous function; i.e., 

f(t,h) = \f{tV* t hW»). (A2) 
b 

This entails the following behaviour of the "magnetisation" 

df(t,h) 



M{t,h) :-- 



(A3) 



dh 

M(t,h) = b Vh - 1 M{tb y \hb yh ) . (A4) 



t fixed 

,-1 



The scaling parameter, b, is arbitrary and along the trajectory \t\b yt = 1 one has 

M(t, h) = It^-^M M (sgn(t), h \t\- yh/yt ) , (A5) 
M(t,0) oc \tf, (3:= (A6) 

yt 

Alternatively, along the trajectory hb Vh = 1 

M(t, h) = h^- yh)/yh M(t hr yttVh , 1) , (A7) 
M(0, h) oc h 1/s , 5 := . (A8) 

This defines the critical behaviour and provides that direct means of extracting the critical 



exponent (3 employed in Refs. jll]jr^] . However, because of numerical noise, it can be difficult 



to extract quantitatively accurate results using this method. 
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The critical exponents can also be determined by studying the pseudocritical behaviour 
of the chiral and thermal susceptibilities, defined respectively as 



Xh{t,h) :-- 
Xt(t,h) :-. 



dM(t, h) 



dh 
dM(t, h) 



Of 



t fixed 



h fixed 



(A9) 
(A10) 



For convenience, we often use Xt{T, h) : = (1/T C ) Xt(t, h). 
For t, h — > + , along hb Vh = 1, one has 

Xh(t,h) = h^^ X h(th- Vt ^,l), 
Xt (t, h) = h^-" t)lyh X t(th- yt/yh , 1) . 



(All) 
(A12) 



At each h, Xh{t,h) and Xt(t,h) are smooth functions of t. Suppose they have maxima 
at £p C and i* respectively, described as the pseudocritical points. Consider the chiral 



susceptibility. At its maximum 



''pc 

d 



dt V 



d_ 

dz 



Xh{z, 1) 



z =th- y t/yh 



(A13) 
(A14) 



f' 1 



which entails that 



tt = K h h yt/yh = Kh h l,m 



pc 



where Kh is an undetermined constant. Similarly, 

t* = K t h yt/yh = K t h 1/m . 



(A15) 



(A16) 



Since (35 > 0, it follows that the pseudocritical points approach the critical point, t = 0, as 
h — > + . It follows from Eqs. (|A15|) and (|A16|) that at the pseudocritical points 



pc 

X h 

pc 

Xt 



Xh (t h pc , h) oc h- Zh 
Xt {tl c ,h)^h- Z \ 



,z h :=l-h 



i 

6 ■ 



(A17) 
(A18) 



Thus by locating the pseudocritical points and plotting the peak-height of the suscepti- 
bilities as a function of h one can obtain values of T c , (3 and 5. 
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TABLES 



logfr Tl l c/v xl c h T pch Xt_ 

-5.0 0.15921 707.0 0.15917 248.5 

-4.3 0.15931 241.9 0.15920 145.4 

-4.0 0.15941 152.9 0.15923 115.3 

-3.3 0.15990 52.19 0.15939 67.33 

-3.0 0.16034 32.91 0.15953 53.34 

-2.3 0.16268 11.31 0.16052 30.91 



TABLE I. The pseudocritical points and peak heights for the chiral and thermal susceptibilities 
in the infrared-dominant model, obtained from Eqs. and ( |35[) respectively. 



IR Dominant UV Improved 0(4) Lattice 

5 3l) 4.3 ± 0.3 4.82 ± 0.05 4.3 ± 0.5 

j3 0.50 0.46 ± 0.04 0.38 ± 0.01 0.30 ± 0.08 

-k 0.67 0.54 ± 0.05 0.55 ± 0.02 0.77 ± 0.14 

T c (MeV) 168.7 153.5 ± 0.1 140 .. . 160 

TABLE II. Critical exponents and temperature for the models considered herein and a compar- 
ison with the results in the N = 4 Heisenberg magnet Q, labelled as 0(4), and lattice simulations 
of QCD with two light flavours ||. 



log h 


T£ c (GeV) 


XT (GeV) 


(GeV) 


pc 

Xt 


-4.30 


0.15464 


896.3 


0.15379 


67.72 


-4.00 


0.15515 


530.7 


0.15394 


55.70 


-3.70 


0.15571 


303.8 


0.15422 


45.70 


-3.52 


0.15627 


224.9 


0.15443 


40.65 


-3.40 


0.15677 


181.8 


0.15460 


37.37 


-3.30 


0.15729 


154.9 


0.15487 


35.00 


-3.15 


0.15795 


120.3 


0.15508 


31.64 


-3.04 


0.15840 


97.21 


0.15534 


29.32 


-3.0 


0.15872 


90.03 


0.15536 


28.39 



TABLE III. The pseudocritical points and peak heights for the chiral and thermal susceptibil- 
ities in the ultraviolet-improved model, obtained from the numerical solution of Eq. (||) with the 
gluon propagator of Eqs. ([2^)-(p4]). From the dependence of the peak heights on the number of 
points in the p-array, we estimate a systematic 1.5% error in and 10% in Xt- 
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FIGURES 
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FIG. 1. The chiral susceptibility, Xh(T,h), in the infrared-dominant model, Eq. (|34|), as a 
function of T for various values of h. 
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FIG. 2. The thermal susceptibility, xt{T,K), in the infrared-dominant model, Eq. (35), as a 
function of T for various values of h. 
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FIG. 3. The chiral susceptibility, Xh(T,h), in the ultraviolet-improved model, Sec. IIB , as a 
function of T for various values of h. 
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FIG. 5. The peak heights at the pseudo critical points of the chiral and thermal susceptibilities 
in the ultraviolet-improved model: X/T (filled-circles), x^r (filled squares). The solid lines are 
straight-line fits, with the slopes 
in Eqs. ([AT7D and ( flip . 



and 



given in Eq. (|39|), which verify the scaling laws 
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